Bulk viscosity of spin-one color superconductors with two quark flavors 
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We consider the contribution of the Urca-type processes to the bulk viscosity of several spin-one 
color-superconducting phases of dense two-flavor quark matter. In the so-called transverse phases 
which are suggested to be energetically favorable at asymptotic densities, the presence of ungapped 
quasiparticle modes prevents that spin-one color superconductivity has a large effect on the bulk 
viscosity. When all modes are gapped, as for one particular color-spin-locked phase, the effect on 
the viscosity can be quite large, which may have important phenomenological implications. 

PACS numbers: 12.38.-t, 12.38.Aw, 12.38.Mh, 26.60.+C 



I. INTRODUCTION 

Dissipative processes play an important role in the evo- 
lution of neutron stars. These processes are governed by 
transport coefficients, such as the heat and electrical con- 
ductivities, the neutrino diffusion coefficient, as well as 
the shear and bulk viscosities. The conductivities are 
important for stellar cooling as well as for the magnetic 
field decay. Shear viscosity dampens differential rotation 
in a star and, thus, leads to a uniform rigid-body rota- 
tion. Bulk viscosity, on the other hand, dampens density 
oscillations inside the star. Both differential rotation and 
oscillations could be excited in newly formed (hot) neu- 
tron stars, or could develop in old (cold) stars due to ex- 
ternal perturbations, e.g., such as matter accreted from 
a companion star. 

It is interesting to note that, in the absence of viscosity, 
all rotating stars would be unstable. The reason is that 
such stars spontaneously develop instabilities as a result 
of the emission of gravitational waves JEilfll (for 
reviews on this topic see, e.g., Refs. 0,01). The so-called 
r-mode (or rotation-dominated) instabilities might be the 
most important ones. They can develop at a relatively 
low angular velocity Q, and therefore may be relevant 
for a large number of compact stars. 

The main theoretical uncertainty in predicting whether 
the r-mode instabilities develop in a star lies in the poor 
understanding of the viscosity of dense baryon matter, as 
well as in the limited knowledge of the stellar composi- 
tion. There is hope, however, that a systematic approach, 
based on a broad understanding of various properties of 
dense baryonic matter, can eventually result in a clear 
picture regarding the neutron star composition. 
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The viscosity of nuclear and mixed phases of dense 
baryonic matter has been calculated under various con- 
ditions and assumptions over the last three decades 
[! M El HI 0, M HE E H El- The bulk viscos- 
ity of normal conducting strange quark matter was also 
calculated The latter might be relevant if the 

baryon density in the central regions of neutron stars is 
so high that matter becomes deconfined. 

The physical conditions in the interior of such stars are 
quite unique: this is the only place in the Universe where 
a deconfined state of cold and dense baryonic matter can 
naturally exist. This possibility has attracted a lot of 
attention since the notion of quarks was introduced pll . 

If deconfined quark matter does exist inside stars, it is 
most likely color-superconducting. (For reviews on color 
superconductivity, see Refs. [H, [27|, S [H HJ SH H2-) 
It is therefore of great interest to study various transport 
properties of color-superconducting phases of quark mat- 
ter. First attempts have already been made to estimate 
the heat and electrical conductivity [H, [34]], as well as 
the bulk and shear viscosities [35S, 36] in the color-flavor- 
locked (CFL) phase of quark matter. Also, in the case 
of the two-flavor color-superconducting (2SC) phase, one 
can argue that most transport coefficients are dominated 
by the two ungapped (blue) quasiparticles (3^. (For a 
recent detailed study of the bulk viscosity in the 2SC 
phase, see Ref. [371]. ) 

In this paper, we calculate the bulk viscosity of the 
four most popular spin-one color-superconducting phases 
of two-flavor (non-strange) quark matter: the color-spin- 
locked (CSL) planar, polar, and the A phase [38l l39l, liol. 
SH, S3, El E3- One of these is likely to be the ground 
state of dense baryon matter if the spin-zero Cooper pair- 
ing of quarks is prevented by the constraints of charge 
neutrality and /3-equilibrium. Moreover, cooling calcula- 
tions for neutron stars favor small gaps of the order of 
1 MeV [45[ which is the typical size of the gap in spin- 
one color superconductors [Hi, [II, S3, Sj], El EE SI ■ The 
absence of strange quarks in the system may be natural 
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if the medium-modified constituent value of the strange 
quark mass is larger than the corresponding value of the 
chemical potential. The generalization of this study to 
the case of spin-one color-superconducting strange quark 
matter will be reported elsewhere 46 1. 

As is well known, in fully gapped spin-zero color- 
superconducting phases, the thermal densities of the 
quark- and hole-type quasiparticles are suppressed ex- 
ponentially by the energy gap, n qp oc exp(— 4>a/T) where 
4>o and T are the values of the spin-zero gap and the tem- 
perature, respectively. This then translates into an ex- 
ponential suppression of the quasiparticle contributions 
to the transport coefficients. In Ref. 35] such an argu- 
ment was used in order to get a simple estimate of the 
viscosity in the CFL phase. It should be noted, how- 
ever, that many transport properties are not dominated 
by the quasiparticles when there exist Nambu-Goldstonc 
excitations in the low-energy spectrum, as is the case in 
the CFL phase [H, [34|, [3(| . The situation is expected to 
be different also in the case of spin-one phases which, in 
general, are not isotropic, and whose gap functions may 
have nodes for some directions of momenta. 

The effect of non-isotropic gaps and various topologies 
of nodes on the neutrino emission and the cooling rate 
of spin-one color-superconducting phases were recently 
discussed in detail [43, HgJ . Following a similar approach, 
in this work we study the bulk viscosity. 

Since the order parameter in spin-one color supercon- 
ductors breaks rotational invariance, dissipative hydro- 
dynamics is more complicated than for isotropic media. 
For instance, we expect that the viscosity becomes a ten- 
sor [49(. We shall avoid these complications by making 
the implicit assumption that the hydrodynamic equa- 
tions are averaged over solid angle. In this way, only 
one (angular-averaged) bulk viscosity coefficient, will 
appear in the hydrodynamic equations and will have to 
be extracted from the (angular-averaged) neutrino emis- 
sion rate. 

In addition, one of the spin-one superconducting 
phases studied here, namely the CSL phase, breaks 
baryon number, i.e., it is also a superfluid. Fortunately, 
it is not an anisotropic superfluid because the order pa- 
rameter does not break rotational invariance. (Here we 
assume that the magnetic field of the star is not strong 
enough to align the spins of the Cooper pairs.) Neverthe- 
less, isotropic superfluids still have a rather complicated 
hydrodynamic behavior, invol ving three (instead of only 
one) bulk viscosity coefficients [50Ll5l|. However, it is not 
completely unrealistic to assume that the relative veloc- 
ity between normal and superfluid components is neg- 
ligible compared to the absolute velocity of the normal 
component. In this case, only one coefficient contributes 
to energy dissipation. In this sense, our treatment is 
completely analogous to that of Ref. [l6| , the difference 
being that here we consider quark matter instead of nu- 
cleonic matter. Let us finally note that the dissipative 
hydrodynamics of anisotropic superfluids is even more 
complicated, see for instance the case of superfluid He-3 
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FIG. 1: Diagrammatic representation of the weak processes 
that contribute to the bulk viscosity of non-strange quark 
matter in stellar cores. 



In the CSL phase the breaking of baryon number gives 
rise to a phonon as the corresponding Goldstone excita- 
tion. We neglect the contribution of the Goldstone mode 
to the bulk viscosity coefficient for the following reason. 
First, the effective theory for phonons is approximately 
scale-invariant at very low energies. For scale-invariant 
superfluids, however, two of the three bulk viscosity co- 
efficients have been shown to vanish [5l[ . The remaining 
coefficient may be non-zero, but corresponds to dissipa- 
tion due to relative motion of superfluid and normal com- 
ponent, which we have already assumed to vanish. 

The remainder of this paper is organized as follows. 
In the next section, we introduce the formalism for cal- 
culating the bulk viscosity in non-strange quark matter. 
In Sec. IIIII we calculate the bulk viscosity in the normal 
phase of two-flavor quark matter. Then, in Sec. IIVI we 
present our results for the bulk viscosity in spin-one color- 
superconducting phases. The discussion of the results is 
given in Sec. fVl 



II. BULK VISCOSITY 

As mentioned in the introduction, the bulk viscosity 
is responsible for the damping of density oscillations in 
compact stars. The characteristic frequencies of interest 
(e.g., set by the r-mode instabilities) are comparable to 
the rotational frequencies of stars, i.e., uj < 10 4 s -1 
These frequencies are many orders of magnitude smaller 
than the typical rates of strong interactions, and there- 
fore quark matter cannot be driven substantially out of 
equilibrium with respect to strong processes. This is the 
reason why the bulk viscosity is dominated by the much 
slower, flavor-changing weak processes [H H3]. In the 
case of non-strange quark matter studied here, the rel- 
evant processes are electron capture by u quarks and (3 
decay of d quarks, see Fig. [TJ 

Let us assume that small oscillations of the quark mat- 
ter density are described by Sn = 5n Re(e lwt ) where 5no 
is the magnitude of the density variations. For such a 
periodic process, the bulk viscosity C is defined as the 
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coefficient in the expression for the energy-density dissi- 
pation averaged over one period, r = 2tt/uj, 



(4iss) = -- / dt(W-vf 



(1) 



where v is the hydrodynamic velocity associated with 
the density oscillations. By making use of the continuity 
equation, h + n V • v = 0, we derive 



diss/ 



C^ 2 / Sn 
2 V n 



(2) 



In order to solve for £, the dissipated energy on the left- 
hand side has to be calculated explicitly. This can be 
done as follows. 

The density oscillations drive quark matter slightly 
out of (3 equilibrium, but not out of thermal equili- 
brium which is restored almost without delay by strong 
processes. The corresponding instantaneous quasi- 
equilibrium state can be unambiguously characterized by 
the total baryon number density n and the lepton frac- 
tion X P , 



1 / 

- (n u + n d ) 

n ' 



(3a) 
(3b) 



where n u and nd are the number densities of up and 
down quarks, while n e is the number density of electrons. 
(In the case of strange quark matter, one should also add 
the strangeness fraction X s = nJn where n s is the num- 
ber density of strange quarks [4y] .) Charge neutrality 
requires 



2 1 
-n u - -n d 



n R = 0. 



(4) 



Using this constraint together with the definitions in 
Eq. ([3]), one can express the number densities and, in 
fact, all thermodynamic quantities of quark matter in 
terms of n and X e . For the number densities, for exam- 
ple, one finds 



n e 
n u 
rid 



X e n, 

{1+X e )n, 
(2 - X e )n. 



(5a) 
(5b) 
(5c) 



These number densities can also be expressed in terms of 
the corresponding chemical potentials, fa = fa(fa). In /? 
equilibrium, the three chemical potentials are related as 
follows: fid = fx u + /i e . In pulsating matter, on the other 
hand, the instantaneous departure from equilibrium is 
described by the small parameter 



Su = /J, d - jU„ - U e = 5u d - 6 flu - Sfi e , 



(6) 



where 8 fa denotes the deviation of chemical potential fa 
from its value in j3 equilibrium. The quantity Sfj, can be 



conveniently expressed in terms of the variations of the 
two independent variables Sn and 5X ei 



5 ii = C— + B5X e 
n 



(7) 



where, as follows from the definition, the coefficient func- 
tions C and B are given by 

dfid d[i u dfa 

C = n d - n u - n e - — , (8a) 

and on u dn e 

B = - n ( — + — + —\ (8b) 
\dn d dn u dn e J ' 

When Sfi is non-zero the two Urea processes, shown 
diagrammatically in Fig. [U have slightly different rates. 
To leading order in 8u, we could write 



-X6u. 



(9) 



(Note that our A is defined so that it is non-negative.) 
The net effect of having different rates for the two pro- 
cesses is a change of the electron fraction in the system: 



d(5X e ) 
1 dt 



XSfa 



(10) 



This has the tendency to restore the equilibrium value of 
X e . Since the rate is finite, however, the weak processes 
always lag behind the density oscillations. In order to see 
this explicitly, we substitute S/i from Eq. (0 into Eq. iflT 
and get the equation for SX e in a closed form, 



d(8X e ) 



dt \ n 



BSX, 



(11) 



The periodic solution to this equation can be found most 
easily by making use of complex variables. Denoting 



SX e = Re (SX e 



o e 



we derive the following result: 



SX efi = 



Sn C 



B 



(12) 



where, by definition, a = nto/X. In the last equation, 
the lagging of the weak processes is indicated by a non- 
vanishing imaginary part of SX e ,Q. Such an imaginary 
part controls the phase shift of the 5X e oscillations with 
respect to the oscillations of density. 

As we show in a moment, the same phase shift also 
leads to a non- vanishing dissipation of the energy density, 



(4 



PVdt 



(13) 



where V = 1/n is the specific volume. 

The pressure oscillations around the equilibrium value 
are driven by the oscillations of its two independent vari- 
ables, i.e., the quark number density and the lepton frac- 
tion, 



dP 

5P = —5n-nCSX e 
on 



(14) 
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where C is the same as in Eq. (I8al) . In the derivation we 
took into account that rij = dP/dpi and that the total 
pressure is given by the sum of the partial contributions 
of the quarks and electrons, P = J^. Pi(fii). 

After taking into account the relation (fT4|) together 
with the solution for 8X e $ in Eq. (Tl2|) . the expression 
(fl~5|) becomes 



we derive 



C 

(£di SS ) = — u>5n lm.(5X ej0 ) 



1 / (Jnr 



Xw 2 C 2 



uS 1 + (XB/nY 



(15) 



By comparing this with the definition in Eq. ([2]), we fi- 
nally derive an explicit expression for the bulk viscosity, 



c 



AC 



u 2 + {XB/nY 



(16) 



This expression shows that the viscosity is maximum in 
the limit of zero frequency, C max = Cw=0i an d that it falls 
off as 1/uj 2 at high frequencies, u> ^S> u>q = XB/n. It 
should be also noted that ojo ~ A, and that the maxi- 
mum viscosity is inversely proportional to the rate, i.e., 
(max ~ 1/A. Because the rates of the weak processes in 
a dense medium usually have a power-law (or even an 
exponential) dependence on the temperature, the bulk 
viscosity is a very sensitive function of the temperature, 
too. 



III. BULK VISCOSITY IN THE NORMAL 
PHASE 

In order to calculate the bulk viscosity in the normal 
phase of two-flavor quark matter, we need to determine 
the corresponding thermodynamic coefficients B and C 
[see Eq. ©] and calculate the difference of the rates of 
the two Urea processes shown in Fig. [T] 

By making use of the following relations valid for non- 
interacting quark matter 



1 



n u .d 



= 4( 



2 



,3/2 



1 

3^ 



(17a) 
(17b) 



C 
B 



3/i u 3u d ' 

,2/1 1 



(18a) 
(18b) 



Here we made use of the equilibrium relation satisfied 
by the chemical potentials, p, d — u u + fi e , and neglected 
higher-order mass corrections in the expression for B. In 
the temperature regime of interest, T < m Ut d, we also do 
not need to take into account any corrections due to a 
non-zero temperature. 

It should be noted that the coefficient C, and there- 
fore the bulk viscosity which is proportional to C 2 , van- 
ishes in the case of massless quarks. Moreover, this state- 
ment remains true even if the following (non-)Fermi liq- 
uid correction due to strong forces are taken into account 

IsilsllsEIES, 



c' = 



3tt 



hi 



2ud 

m d 



/'-« 



In 



2u v 



(For a recent discussion of (non-)Fcrmi liquid corrections 
see, for example, Ref. [5q.) Because of the large loga- 
rithms on the right-hand side of Eq. (fT9|) . this correction 
can become even larger than the leading-order result for 
C in Eq. (|18ap . Our estimates show that, in two-flavor 
quark matter, taking C into account may increase the 
value of the viscosity by approximately an order of mag- 
nitude. Therefore, in all numerical estimates below we 
add the contributions from Eqs. (|18a[) and (fn)|) . 

Now let us turn to the calculation of A define d by 
Eq. © . Following the original approach of Iwamoto [53] , 
we get the rate for (3 decay in the following form: 



T (6n) = 6 



d 3 p d d 3 p u d 3 p p d 3 p e 
{2-KfE d E u E D E e 



\M\ 2 6 i (P d -P u -P e - P 9 ) f (E d - W ) [1 - / {E u - Mu)] [1 - / (E e - Me)] ■ (20) 



r 



Here, Pi and pi are the 4- and 3- momenta of the «th squared is given by 57] 
particle, respectively, and f(E) = \/{e E / T + 1) is the 

Fermi distribution function. The scattering amplitude \M\ = 64G F cos 9c(Pd ■ Pv){Pu • Pe) 

w -7T -G 2 F cos 2 9c E u E d E P E e (l - cos^ p )(21) 

37T 
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FIG. 2: The bulk viscosity for the normal phase of two-flavor 
quark matter as a function of the period of the density oscil- 
lations. 



After substituting this approximate form of \M\ 2 , all an- 
gular integrals in Eq. ([20)) can be done exactly. Then, 
using the dimensionless variables Xi — [Ei — fii)/T (note 
that [ly = 0), we obtain 

4a s r f°° 

= —&-Gp cos 2 OcHd^uHeT^ I dx v x\J{xv - £), 
n Jo 

(22) 

where £ = dfi/T, and 



J(x) 



n 

;=i 

^2 



dxjf(xj) 



S(xi + x 2 + x 3 - x) 



2(l + e x )' 



By noting that 
expression 



(23) 

r K (— £), we finally derive the 



A = 



17 



G 2 F cos 2 8ca s fj,d^u^eT 4 



17. 
15 



G 2 , cos 2 eca^Xey^nT 4 [1 + 0(X e )] .(24) 



This result, together with the explicit form for the co- 
efficient functions C, C", and B in Eqs. (fT8|) and (fl"9|) . 
is sufficient to calculate the bulk viscosity of the normal 
phase of dense quark matter. The numerical results are 
presented in Fig. [51 In the calculation, we used the fol- 
lowing representative values of the parameters: 

fi d = 500 MeV, m d = 9 MeV, (25a) 
fi u = 400 MeV, m u = 5 MeV, (25b) 
/i e = 100 MeV, a s = 1. (25c) 



In order to understand the numerical results better, 
we calculate the maximum value of the bulk viscos- 
ity, Cmax = C( w = 0), an d the characteristic frequency 
ma = XB /n, which separates the low- and high-frequency 
regimes, 



7.26 x 10 2 



T V 



" * 4,) "fiw) s 1 



— J gcm-s-, (26) 
(27) 



The knowledge of these two quantities is sufficient for 
determining the bulk viscosity at all frequencies, 



Cn 



1 + (lu/lu )< 



(28) 



Taking into account the relations in Eqs. (|26| and (|27p . 
the temperature dependence of the results in Fig. [2] be- 
comes clear. In accordance with the first relation, the 
maximum value of the viscosity (approached at low fre- 
quencies or large r) becomes larger with decreasing tem- 
perature. Also, in agreement with Eq. (|2"?| . the location 
of the "shoulder" (occurring at To = 2tt/ojq) shifts to 
smaller periods (higher frequencies) when the tempera- 
ture gets higher. 



IV. BULK VISCOSITY IN SPIN-ONE 
COLOR-SUPERCONDUCTING PHASES 



In this section, we calculate the bulk viscosity of several 
spin-one color-superconducting phases that have been 
proposed and studied in detail in Refs. [39. |4[| |42|. |43|. |44|| . 
Following Ref. [47] . we focus on the so-called transverse 
phases in which only quarks of opposite chiralities pair. 
Theoretical studies at asymptotically large densities sug- 
gest that such phases are preferred [4fJ, [4J] . 

In (3 equilibrium, the rates of the Urea processes were 
calculated in Ref. (47| in four different spin-one phases. 
Following the same method, here we generalize the cal- 
culation to the quasi-equilibrium state, characterized by 
a small but non- vanishing value of Sfi = fid — fi u — fi e . 

Starting from equations analogous to Eq. (36) and (37) 
in Ref. [47]], we derive 
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(29) 



where S = cos# u = cos Od is the angle between the three-momentum of the u/d quark and the z-axis. The functions 
F^ Vd (E, x) are the same as in Ref. 47J. Their explicit form is given by 



F£ {R,x) = w„(E) V / / dx d dx u (e-^ 

e 1 ,e 2 =± J V J V 
X (^ e x "+ e ^^/ x l J r X =,rVl-e-2^Jx' 2 d +\^^ip 2 d _|_ j 



(30) 



where, by definition, £ = 5fx/T and ^ = 4>i/T. For an 
explicit form of the functions As ir . and w rr , we refer the 
reader to Refs. [ii.!^. 

By expanding the rate difference (|29p in powers of £ 
and extracting the coefficient of the linear term, we derive 
an expression for A in the following general form: 



\{<Pd,<fiu) = A(0) 



12, 

3 + ^H{tp d ,ip u ) 



(31) 



where A(0) is the same as in the normal phase of quark 
matter, see Eq. (|24|) . and H(ipd, <p u ) is a reduction factor 
whose explicit form depends on the choice of a specific 
spin-one color-superconducting phase. Note that, as in 
the case of the neutrino-emission rates [U Hsj , A consists 
of two qualitatively different contributions. The first one 
is given by the term 1 /3 in the square brackets of Eq. (f3Tj) . 
It originates from the ungapped modes that are present in 
all considered spin-one phases. The second contribution 
is given by the term proportional to H((p u , tfd)- This one 
originates from the gapped modes. An explicit form of 
the function H{Lp u ,ipd) in the case of the CSL, planar, 
and polar phases is given by 



60 
17tt 4 



dE 

-1 JO 



ii 



and, in the case of the A phase, it reads 



(32) 



60 
17tt 4 



dE 
i Jo 



dx,. 



d (E,x u )-F^ ulpd (E,x u )]. (33) 



In the derivation we integrated by parts so that the re- 
sults are given in terms of the function (S, x). The 
reduction factors can be easily calculated numerically for 
each phase. The results are shown in Fig. [3] for the case 
fu = fd = V- In the limit of large <p analytical expres- 
sions for the suppression functions for all four spin-one 
color-superconducting phases are given in Appendix O 

For many applications, it is of interest to know the tem- 
perature dependence of the rates. Therefore, in Fig. 21 
we also show the temperature dependence of the suppres- 
sion factor H(T). In the calculation we used the following 



3. 
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FIG. 3: The reduction factor as a function of = <j>/T for 
the CSL, planar, polar, and A phases. 



temperature dependence of the gap parameter, 



<KT) = <t>o \ 1 



(34) 



with (f>o being the value of the gap parameter at T = 0, 
and T c being the value of the critical temperature. Note 
that the ratio T c / <fio depends on the choice of the phase 
The approximate values of this ratio are 0.8 (CSL), 
0.66 (planar), 0.49 (polar), and 0.81 (A phase). 

With the expression for A at hand, we are now in the 
position to calculate the bulk viscosity of the spin-one 
color-superconducting phases. Before we do this, it might 
be appropriate to mention that the coefficient functions 
B and C do not change much due to superconductiv- 
ity. Corrections to B are of order {4>i / Hi) 2 and, thus, 
are strongly suppressed. The function C, on the other 
hand, gets corrections of order (fq / fii which are negligi- 
ble only if <j)f <C m 2 . We assume that this is indeed the 
case. Notably, the spin-one gap corrections to C should 
be comparable to the corrections due to a non-zero tem- 
perature that we neglected in our calculations. 

It should be emphasized that the bulk viscosity in the 
color-superconducting phases has the same general struc- 
ture as in the normal phase, see Eq. (|28|) , but the quanti- 
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FIG. 4: The reduction factor as a function of the temperature 
for the CSL, planar, polar, and A phases. 

ties Cmax and loq should be redefined to take into account 
the rate suppression factors: 

, S pl _ Cmax ,or\ 

Wax — 7 i \°°) 
"spl 

Wg pl = h spl ujQ, (36) 

where C max and ojq are given in Eqs. (f2"6")l and Q27p. re- 
spectively, and 

h sp i = \ + \h(T/T c ). (37) 

We see that the effect of the suppression of A due to 
Cooper pairing manifests itself in a non-trivial way in 
the expression for the bulk viscosity, i.e., 

j. _ Cmax^spl Cro\ 

Qpl - h* pl + (to/corf- W 

From the analysis of this representation, we find that the 
suppression of the rates tends to decrease the viscosity at 
high frequencies (u > /ispi^o) and to increase it at low 
frequencies (oj < h sp iLu ). One should note, though, that 
the relevant range of low frequencies would shrink a lot 
if the suppression happened to be strong. 

The representation in Eq. (|3S[) shows that the effect 
of color superconductivity cannot be very large. Even if 
the suppression of the weak rates due to gapped modes is 
maximal, i.e., H(T/T C ) = 0, the results for Cmax and ojq 
could not change by more than a factor of 3 compared 
to the normal phase of matter. Of course, this is the 
consequence of having ungapped quasiparticle modes in 
the energy spectra. 

The numerical results for the bulk viscosity in the 
normal phase and spin-one color-superconducting phases 
with different values of the critical temperature are shown 
in Fig. [5] For each of the three choices of T c , we show 
a shaded area (color online: blue, green and red for 
T c = 0.25, 1, and 4 MeV, respectively), that is bounded 
by the values of the bulk viscosity in the CSL phase (solid 



FIG. 5: The bulk viscosity as a function of the temperature for 
spin-one color superconducting quark matter. The oscillation 
frequency is u) = 10 3 s _1 . 




0.1 0.2 0.5 1 2 5 



T [MeV] 

FIG. 6: The bulk viscosity as a function of the temperature 
for a toy model of a spin-one color superconductor in which 
all quasiparticle modes are gapped. The oscillation frequency 
is u = 10 3 s _1 . 



boundary) and in the A phase (dotted boundary). In 
the calculation, we used the set of parameters given in 
Eq. (25). 

For comparison, in Fig. [6] we also show the results for 
the bulk viscosity in a toy model in which all quasiparti- 
cles modes are gapped. In this case, the effect is indeed 
very dramatic. The largest deviations from the result of 
the normal phase correspond to the CSL phase (thin solid 
lines), and the smallest to the A phase (dotted lines). 
The viscosity in the planar phase (long-dashed lines) is 
numerically close to that in the CSL phase, and the vis- 
cosity in the polar phase (short-dashed lines) is typically 
somewhere in between. 

At this point it is natural to ask ourselves if there exists 
the possibility that all quasiparticles in a spin-one color 
superconductor are gapped. The answer is affirmative. 
This is the case, ejr., in the version of the CSL phase 
proposed in Ref. [58j |. when the quarks are massive. In 
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our notation, the three quasiparticle gaps are given by 



T>i,2 



,3 - —f=—, (i = u,d). (39) 



v 7 ^ 



The qualitative change in the low-energy spectrum comes 
from a different choice of the gap matrix in the CSL 
phase. (For a discussion of the specific differences, see 
Sec. VII in Ref. [47| ■) In fact, it is likely that the ansatz 
that produces the fully gapped phase is energetically fa- 
vored [5!| • The most important effect of the additional 
non-zero gap is the appearance of another suppression 
factor multiplying the first term in the expression for A, 
see Eq. (|3~Tj) . Therefore, the modified expression for the 
CSL phase reads 



X(ipd, <p u ) = A(0) 



1 „ / m d (p d m u ip u 
3 Vv^Md' V2^„ 



;H((p d ,ip u ) 



(40) 

Note that the additional suppression factor is given in 
terms of the same function H(ip,i,y>u) which was cal- 
culated numerically for the transverse CSL phase, see 
Figs. [3] and 2J If the up and down quark masses are 
rather small, the effect from the additional suppression 
cannot be easily seen before the temperature becomes 
much lower than the critical value, i.e., T -C T c . The 
situation changes dramatically, however, if the relevant 
constituent values of the quark masses happen to be con- 
siderably larger than the current masses of quarks. As 
suggested by the analysis in Ref. [fill , this is indeed pos- 
sible. Then, the bulk viscosity could be affected almost 
as much as in the toy model in Fig. [6] 



V. DISCUSSION 

In this paper we have calculated the bulk viscosity 
for the normal phase as well as for four spin-one color- 
superconducting phases of two-flavor dense quark mat- 
ter. The main contributions come from the Urea pro- 
cesses shown diagrammatically in Fig. [TJ Note that the 
results for the normal phase are also relevant for the 2SC 
phase. Indeed, after taking into account that there are 
two (blue) ungapped modes of quasiparticles in the low- 
energy spectrum of the 2SC phase, the low-temperature 
bulk viscosity is approximately given by the same expres- 
sion H28D . provided the following redefinitions are made: 



f 2SC _ n C 



and ujq = u>o/3, where the normal-phase 



quantities are given in Eqs. ([26]) and (f27j) . respectively. 
The redefinitions account for the decrease of the weak 
rates by a factor of 3 at T <C Ao where Ao is the value 
of the 2SC gap. 



The microscopic calculations of the bulk viscosity in 
the spin-one color-superconducting phases suggests that 
quasiparticles with different types of gapless nodes (e.g., 
points or lines at the Fermi sphere) could potentially 
play a very important role. In the case of the trans- 
verse phases, however, the presence of a single ungapped 
quasiparticle mode washes out essentially all information 
about spin-one Cooper pairing, see Fig. [fTJ The presence 
of non-zero quark masses may provide a gap for such a 
mode and the situation changes. In this paper, we briefly 
discussed such a possibility in connection with the CSL 
phase of Ref. [58[. The results are shown in Fig. [51 

In agreement with the general expectation, we find 
that the bulk viscosity often tends to decrease when there 
is Cooper pairing of quarks whose main effect is to sup- 
press the rates of the weak processes. In some cases (e.g., 
at sufficiently low frequencies and/or at temperatures 
close to the critical value) the behavior may reverse be- 
cause of the non-trivial dependence of the bulk viscosity 
on the suppression factor, see Eq. (|38[> . Such an increase 
of the viscosity in the color-superconducting CSL phase 
is seen, for example, in a range of temperatures below T c 
in Figs. [5] and [6] in the case when T c = 4 MeV. 
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APPENDIX A: ASYMPTOTIC FORM OF THE 
SUPPRESSION FUNCTION H(tp,<p) AT ip -> oo 

In this appendix, we calculate the asymptotic form of 
the suppression function H(tp, ip) at ip — > oo. This is 
relevant for the bulk viscosity of spin-one color super- 
conductors in the low-temperature regime. 

Let us start from the CSL and planar phases whose 
gap functions have no nodes in momentum space. To 
this end, it is useful to compute the asymptotic behavior 
of the following integral: 
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OC /> OO /» OC 



J{(p,<p) = V / / / dx u dx d dx,xJe-^V^^ + l) 'L^V^l^ + l 

e u e 2 = ± J JO JO V ' V > 



By neglecting the terms of order 0(e 2v ), we arrive at 
the following leading-order asymptotic behavior: 



J{ip,(p) ~ 4<^e v 

/o "'O 
5.047(^e" v 



oc roc />oo 



dy u dy d dx v x v 



(A2) 



Note that in the derivation we changed the integration 
variables, Xi = yi^/2ip for i = u,d, and calculated the 
remaining integral numerically. 

By making use of the result for J(<p, ip), we derive the 
large tp asymptotic behavior for the suppression function 
H(ip, ip) in the CSL phase: 



240 
17tt 4 



J(V2ip, V2ip) ~ 1.034^-^. 



(A3) 

Similarly after taking into account the angular depen- 
dence of the gap in the planar phase, we derive 

iP lanar (^) = y 1 dsj(Vi + sWi + sv) 

~ 0.917^e^. (A4) 



r 



The derivation in the polar and ^4 phases is slightly 
more complicated because the corresponding gap func- 
tions have zeros for some directions in momentum space. 
By approximating the angular integrals in the same way 
as in Ref. [47| (see Appendix E there), we arrive at the 
following asymptotic behaviors: 



# polar (<P,<p) 
H A (<p,v) 



1 



(A5) 
(A6) 



In connection to these last two results, we should men- 
tion that while the parametric dependence on <p is easy 
to extract analytically, it is much harder to determine 
the overall coefficients in their power-law asymptotic be- 
havior. In our derivation, therefore, we combined the 
analytical derivation with the numerical calculations. 
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